The Modulational Instability for a Generalized Korteveg-DeVries 

equation 



We study the spectral stability of a family of periodic standing wave solutions to the generalized 
KdV (g-KdV) in a neighborhood of the origin in the spectral plane using what amounts to a rigorous 
Whitham modulation theory calculation. In particular we are interested in understanding the role 
played by the null directions of the linearized operator in the stability of the traveling wave to 
perturbations of long wavelength. 

A study of the normal form of the characteristic polynomial of the monodromy map (the periodic 
Evan's function) in a neighborhood of the origin in the spectral plane leads to two different instability 
indices. The first index counts modulo 2 the total number of periodic eigenvalues on the real axis. 
This index is a generalization of the one which governs the stability of the solitary wave. The second 
index provides a necessary and sufficient condition for the existence of a long-wavelength instability. 
This index is essentially the quantity calculated by Haragus and Kapitula in the small amplitude limit. 
Both of these quantities can be expressed in terms of the map between the constants of integration 
for the ordinary differential equation defining the traveling waves and the conserved quantities of the 
partial differential equation. These two indices together provide a good deal of information about 
the spectrum of the linearized operator. 

We sketch the connection of this calculation to a study of the linearized operator - in particular 
we perform a perturbation calculation in terms of the Floquet parameter. This suggests a geometric 
interpretation attached to the vanishing of the modulational instability index previously mentioned. 
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1 Introduction and Preliminaries 

In this paper, we consider standing wave solutions to the generalized KdV (gKdV) equation 

u t = u xxx + (f(u)) x - cu x (1) 



where /(•) G C 2 (M) is a prescribed nonlincarity and c is the wavespeed. Such solutions represent traveling 
wave solutions to the generalized KdV equation with nonlincarity f{u). Of particular interest is the case 
of power law nonlinearity f(u) = u p+1 , which in the cases of p =1,2 represents the equations for traveling 
wave solutions to the KdV and MKdV, respectively. Obviously such traveling wave solutions are reducible 
to quadrature: they satisfy 

(2) 

au + E. (3) 

We are interested in the spectrum of the linearized operator (in the moving coordinate frame) 

H v = v xxx + {vf'(u)) x - cv x 

in two related settings. First, we study the spectrum in a neighborhood of A = 0. Physically this 
amounts to long-wavelength perturbations of the underlying wave profile: in essence slow modulations of 
the traveling wave. There is a well developed physical theory, commonly known as Whitham modulation 
theory [40, 41], for dealing with such problems. On a mathematical level the origin in the spectral plane 
is distinguished by the fact that the ordinary differential equation giving the traveling wave profile is 
completely integrable. Thus the tangent space to the manifold of traveling wave profiles can be explicitly 
computed, and the null-space to the linearized operator can be built up from elements of this tangent 
space. We show that these considerations give a rigorous normal form for the spectrum of the linearized 
operator in the vicinity of the origin providing that certain genericity conditions are met. Assuming that 
these genericity conditions are met we are able to show the following: there is a discriminant A which can 
be calculated explicitly. If this discriminant is positive then the spectrum in a neighborhood of the origin 
consists of the imaginary axis 1 with multiplicity three. If this discriminant is negative the spectrum of 
the linearization in the neighborhood of the origin consists of the imaginary axis (with multiplity one) 
together with two curves which leave the origin along lines in the complex plane, implying instability. 
Long wavelength theories are invariably geometric in nature, and the one presented here is no exception: 
both the instability index and the genericity conditions admit a natural geometric interpretation. 

Secondly, we are interested in determining sufficient conditions for the existence of unstable spectrum 
supported away from A = 0. Here, this is accomplished by calculating an orientation index using Evans 
function techniques: essentially comparing the behavior of the Evans function near the origin with the 

1 Note that this does not imply spectral stability since there is the possibility of bands of spectrum off of the imaginary 
axis away from the origin. 



Uxx + f(u) -CU = 
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asymptotic behavior near infinity. Physically, such an instability amounts to an instability with respect to 
finite wavelength perturbations. The derived stability index is a generalization of the one which governs 
stability of solitary waves. In fact, in the case of power-law nonlinearity and wave speed c > 0, we show 
that in a long wavelength limit the sign of this index, which is actually what determines stability, agrees 
with the sign of the solitary wave stability index derived by, for example, Pego and Weinstein[35, 34]. 

This paper uses ideas from both stability theory and modulation theory, and thus there is an ex- 
tensive background literature. Most obviously is the stability theory of solitary wave solutions to KdV 
type equations which was pioneered by Benjamin[4] and further developed by Bona[6], Grillakis[18], Gril- 
lakis,Shatah and Strauss [19], Bona, Souganides and Strauss [7], Pego and Weinstcin[35, 34], Weinstein[38, 
39] and others. In this theory the role of the discriminant is played by the derivative of the momentum 
with respect to wave-speed. Our discriminant is somewhat more complicated, which is to be expected: the 
solitary waves homoclinic to the origin are a codimension two subset of the family of periodic solutions, 
so one expects that the general stability condition will more complicated. There are also a number of 
calculations of the stability of periodic solutions to perturbations of the same period, or to perturbations 
of twice the period, due to Angulo Pava[l], Angulo Pava, Bona and Scialom[2] and others. In this setting 
the linearized operator has a compact resolvent, so the spectrum is purely discrete, and the arguments 
are similar in spirit to those for the solitary wave stability. In contrast we consider the case of a general 
L2 case where one must understand the continuous spectrum of the operator. 

A stability calculation in the spirit of modulation theory was given by Rowlands [36] for the cubic 
nonlinear Schrodingcr equation. Other stability calculations in the same spirit, but differing greatly in 
details and approach, were given by Gallay and Haragu§[15, 16], Haragus, and Kapitula[21], Bridges and 
Rowlands[ll], and Bridges and Mielke[10]. The work of Gardner[17] is also related, though it should be 
noted that the long-wavelength limit in Gardner is very different from the one we consider here: in the 
former it is the traveling wave itself which has a long period. In our calculation the period is fixed and 
we are considering perturbations of long period. The current paper also owes a debt to the substantial 
literature on Whitham theory for integrable systems developed by Lax and Levermore[28, 26, 27], Flashka, 
Forest and McLaughlin[14], and many others. We note, however, that the calculation outlined in this 
paper is not an integrable calculation. The papers that are perhaps closest to that presented here are those 
by Oh and Zumbrun[31, 32, 33] and Serre[37] on the stability of periodic solutions to viscous conservation 
laws, where similar results relating the behavior of the linearized spectral problem in a neighborhood of 
the origin to a formal theory of slow modulations are proved. 

Our results are most explicit in the case of power law nonlinearity. It should be noted that due to the 
scaling invariance in this case we can always assume that c G {— 1, +1}. Indeed, it is easy to check that 
if u(x,t;c) solves (1) with the nonlinearity f(u) = u p+1 , then 



solves (1) with wave speed sgn(c). 

The paper is organized as follows: in the second section we lay out some basic general properties of 
the spectrum of the linearized operator. In the third section we explicitly compute the monodromy map 
and associated periodic Evan's function at the origin. A perturbation analysis in the neighborhood of 
the origin gives a normal form for the Evan's function. In the fourth section we develop similar results 
from the point of view of the linearized operator: we compute the generalized null-space of the linearized 
operator in terms of the tangent space to the ordinary differential equation defining the traveling wave. 
The structure of this null-space (under some genericity conditions) reflects that of the monodromy map 
at the origin, and a similar perturbation analysis gives a normal form for the spectrum. While the two 
approaches are in principle the same most of the calculations arc more easily carried out in the context 
of the monodromy map/Evans function. We mainly present the calculations at the level of the linearized 
operator since it helps to clarify some aspects of the monodromy calculation. Finally we end with some 
concluding remarks. 




(4) 
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It should be noted we restrict neither the size of the periodic solution nor the period. Moreover, all 
of our analysis applies to both localized and bounded perturbations of the underlying wave. Also in this 
paper "stability" will always means spectral stability. 



1.1 Preliminaries 

Note that the partial differential equation has (in general) three conserved quantities 



M 



= / u(x,t)dx 
Jo 

P= u 2 (x,t)dx 
Jo 

H = J o \ul + F{u)dx 



which correspond to the mass, momentum and Hamiltonian of the solution respectively. These quantities 
considered as functions of the traveling waves parameters will form an important part of the analysis. 

The periodic standing wave solutions of (1) arc of the form u(x,t) = u{x) where u is a periodic 
function in the x-variable. Substituting this into (1) and integrating twice we see that u satisfies 

lu 2 x + F(u)-^u 2 -au = E (5) 

where a, E e R arc constants of integration and F' = f. Note that the solitary wave case corresponds 
to a = 0, E = 0, so the solitary waves are a codimension two subset of the periodic waves. In order to 
assure the existence of a periodic orbits, we must require that the effective potential 

V(u; a, c) = F(u) — ^u 2 — au 

has a local minimum. Note that this places a condition on the allowable parameter regime V for our 
problem. We will always assume that we are in the interior of this open region, and that the roots u+, U- 
of the equation V(u;a, c) = E with V(u;a,c) < E for u e are simple, guaranteeing that they 

are C 1 functions of a, E, c. 

As is standard, the the period of the corresponding periodic orbit is given by 

f u + dii 

T = T(a,E,c):=2 = . (6) 

Ju. y/2 (E-V(u;a,c)) 

The above interval can be regularized at the square root branch points u_, u + by the following proce- 
dure: Write E — V(u; a, c) = (r — u-)(u+ — r)Q(u) and consider the change of variables u — u ++ u - + 
" + ^"~ sin(fl). Notice that Q(u) ^ on [ti_,u + ]. It follows that du = \J(u — U-)(u + — u)d6 and hence 
(6) can be written in a regularized form as 

d6 



T(a,E,c) = 2 J' 



Q(H±+I£_ + H±_i£_ sm (0)) 



Similarly the mass, momentum, and Hamiltonian of the traveling wave are given by the first and second 
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T ' u+ u du 



moments of this density, i.e. 

M(a,E,c) = (u) = I u(x)dx 
Jo 

f T r 

P(a, E, c) = (u 2 ) = / u 2 (x)dx = 2 1 

JO Ju 



^2 (E-V(u;a lC )) 
u 2 du 



_ V2 (E-V(u;a,c)) 

V ' \ 2 V V J u _ y/2(E-V(u;a,c)) 

Notice that these integrals are regularized by the same substitution. In particular one can differentiate 
the above expressions with respect to the parameters (a, E, c) and the derivatives of these quantities will 
play an important role in the subsequent theory. Note that there is a third constant of integration xo 
corresponding to translation invariance, but this can be modded out and does not play an important role 
in the theory. 

These quantities satisfy a number of identities, as is derived in the appendix. In particular if we define 
the classical action 

/r T r u + 
pdq= u 2 x dx = 2 y / 2(E - V(u;a,c)) du 
Jo Ju- 

(which is not itself conserved) then this quantity satisfies the following relations 

K E =T 
K a = M 

K C = P -. 
2 

Using the fact that T, M, P and H are C 1 functions of parameters (a, E; c), the above implies the following 
relationship between the gradients of these quantities 

EVT + aVM + ^VP + VH = 

where V = (<9 a , <9b, 9 C ): see the appendix for details of this calculation. The subsequent theory is 
developed most naturally in terms of the quantities T, M, and P. However, it is possible to restate our 
results in terms of M, P and H using the above identity. This is desirable since these have a natural 
interpretation as conserved quantities of the partial differential equation. 

As noted before this long-wavelength calculation is geometric, and a number of Jacobian determinants 
arise. We adopt the following notation for 2 x 2 Jacobians 



{f,g}x, y = 



fx Sx 

fy 9y 



with {f,g,h} Xi y jZ representing the analogous 3x3 Jacobian. 

We now begin our study of linear stability of the periodic waves u(x) = u(x; a, E, c) under small 
perturbation. To this end, we consider a small perturbation of the periodic wave u(x; a, E) of the form 

ip(x, t; a, E, c) — u(x; a, E, c) + ev(x, t) + 0(e 2 ), 

where < |e| <§; 1 is a small parameter. Substituting this into (1) and collecting the 0(e) terms yields 
the linearized equation d x L[u]v = —Vt, where L[u] := —d 2 — f'{u) + c is a linear differential operator 
with periodic coefficients. Since the linearized equation is autonomous in time, wc may seek separated 
solutions of the form v(x,t) — e^ fJ,t v(x), which yields the eigenvalue problem 

d x C[u]v = [iv. (7) 
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Note that we consider the linearized operator d x C[u] as a closed linear operator acting on a Banach 
space X with domain V(d x C[u\). In literature, several choices for X have been studied, each of which 
corresponding to different classes of admissible perturbations v. In our case, we consider X = L 2 (R;M.) 
and V{d x C[u\) = H 3 (R), corresponding to spatially localized perturbations. In this case standard Floquet 
theory yields the following definitions. 

Definition 1. The monodromy matrix M(/i) is defined to be the period map 

M(/i)=*(T, M ) 

where <&(x, /z) satisfies 

$ x =H(x;/i)$ *(0,/i) = I (8) 

with I the 3 x 3 identity matrix and 

/ 10 

H(x;/i) = 1 

V -M - W(u) -/'(«) + c 

Given the monodromy the spectrum is characterized as follows: 

Definition 2. We say \i e spec(9a;£[u]) if there exists a non-trivial bounded function ip such that 
d x C[u]ip — /j.tjj or, equivalently if there exists a A e C such that |A| = 1 and 

det[M(»-AI] = 0. 

Following Gardner [1 7] we define the periodic Evans function to be 

D(fi, A) = det ( M(/i) — AI) . (9) 

Moreover, we say the periodic solution u(x;a,E,c) is spectrally stable if spec(d x £[u}) does not intersect 
the open right half plane. 

Remark 1. Notice that due to the Hamiltonian nature of the problem, spec(d x C[u]) is symmetric with 
respect to reflections across the real and imaginary axes. Thus, spectral stability occurs if and only if 
sj>ec(d x £[u\) C Mi. Since we are primarily concerned with roots of D(fi,X) with A on the unit circle we 
will frequently work with the function D([i,e lK ), which is actually the function considered by Gardner. 

In this paper, we will study different asymptotics of this function. In the next section, we will study 
the asymptotics of (9) as /x — > oo. This will provide information about the global structure of the 
spectrum of the linearized operator d x C\u\, as well as providing us with a finite wavelength instability 
index which counts modulo 2 the number of intersections of the spectrum with the positive real axis. We 
then study the asymptotics of (9) in the limit (fj,, n) — > (0, 0), which yields a quantity which we refer to as 
a modulational stability index, which is expressed in terms of the derivatives of the monodromy operator 
at the origin. 



2 Global Structure of spec(d x C[u}) 

In this section, we review some basic global features of the spectrum of the linearized operator d x C[u] 
which are useful in a local analysis near \i = 0. We also state some important properties of the Evans 
function D(/j,, A) which are vital to the foregoing analysis. 

Proposition 1. The spectrum spec (d x C[u\) has the following properties: 



2 GLOBAL STRUCTURE OF SPEC(<9 X £[[/]) 



• There are no isolated points of the spectrum. In particular, the spectrum consists of piecewise smooth 
arcs. 

• D(p, A) = det( M(fi) - AI) = -A 3 + a(p)X 2 - a(- f i)\ + 1 with a(n) = tr( M(/x)). 

• The function a(/i) satisfies a(0) = 3, a'(0) = 0. 

• T/ie entire imaginary axis is contained in the spectrum, i.e. iR C spcc(d x £). Further for \fj,\ 
sufficiently large along the imaginary axis the multiplicity is one. 

• IR n spec(d x £) consists of a finite number of points. In particular there are no bands on the real 
axis. 

Proof. The first claim, that the spectrum is never discrete, follows from a basic lemma in the theory of 
several complex variables: namely that, if for fixed A* the function D([i, A*) has a zero of order k at n* 
and is holomorphic in a polydisc about (/i* , A*) then there is some smaller poly disc about (^*, A*) so that 
for every A in a disc about A* the function £>(/i, A) (with A fixed) has k roots in the disc \fi — < S. 
For details see the text of Gunning [20]. It is clear from the implicit function theorem that /j, is a smooth 
function of A as long as ^ = tr(cof*( M(/i) — AI) M^) ^= 0, where cof represents the standard cofactor 
matrix. 

The second claim is an easy symmetry calculation. The stability problem is invariant under the the 
map x — x,[it-^> — (i, which implies that 



M(p) - M _1 (-/x). 



Thus one has 



det[M(-/x)- AI] = det[M _1 (/i)-AI] 

= -A 3 det[M- 1 (/i)]det[M(/i)-A" 1 ] 
= -A 3 (-A -3 + a(/i)A~ 2 + 6(/i)A _1 + l) 
= -A 3 - b(fi)X 2 - a(fi)X + 1 

from which it follows b(fi) = —a(—fx). 

The proof of the third claim will be deferred until lemma 2. 

The fourth claim is another symmetry argument. Since a(/x) is real on the real axis it follows from 
Schwarz reflection that for /i e Mi, we have a(—fi) = a(~p) = a(/z) and the characteristic polynomial takes 
the form 

D(fi, A) = -A 3 + a\ 2 - aX + 1 

where a = a{p), and thus 



— l 



£>(//, A) = -X A D [ijl;X 

Hence for imaginary fi the eigenvalues of the monodromy are symmetric with respect to the unit circle 
with the same multiplicities. Since the monodromy has three eigenvalues, it follows that at least one 
must lie on the unit circle. 

To see that the multiplicity is eventually one we note that by standard asymptotics the monodromy 
M(/i) satisfies 

M(/x) w e A ^ T , |m| » 1 

where A(/x) is defined by 

A(/x) = 
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The three eigenvalues of e A ^ T are given by 

Ai = eS' 3T , A 2 - e-^ T , and A 3 = C -^ T (10) 

where u = e 27 ™/ 3 is the principle third root of unity. If p E R + i it follows that Ai = exp (-|/j| 1 / 3 e" r / 6 T) 
and since cos(7r/6) > we have |Ai| — > as R + i 3 p — > oo. Similarly, A 2 = exp (— |/z| 1 / 3 e 57r / 6 T') and 
A3 = exp (l/il 1 / 3 ?) so that |A 2 | — > 00 as R+i 3 p — > 00 and | A3 1 = 1. Thus, for /U e R + i large, we have 
that p is an eigenvalue of multiplicity one. Similarly, we can show |Ai| — » 00, |A 3 | — > as R+i 9 p — > —00 
and |A 2 | = I for /li € R~i, |/z| ;» 1. Therefore, it follows that /i € spec(9 a; £[M]) with multiplicity one for 
p G R«, ImI > !• 

The final claim follows from a similar asymptotic calculation together with an analyticity argument. 
Notice that for p real the eigenvalues of the monodromy are either all real or one real and one complex 
conjugate pair. If the eigenvalues lie on the unit circle then in the first case 1 or —1 must be an eigenvalue. 
In the second one must have a complex conjugate pair of eigenvalues, and thus (since the determinant is 
one) 1 must be an eigenvalue. Thus if a point on the real axis is in the spectrum then either det( M(/i) — I) 
or det(M(/j) + I) must vanish. Since det(M(/z) ± I) arc entire functions it follows that either they are 
identically zero or the zero set has not finite accumulation points. The large p asymptotics implies that 
they cannot be identically zero, therefore the zero set must be discrete. Further the large p asymptotics 
implies that for sufficiently large p along the real axis p ^ spec (d x £[u]), so the spectrum is confined to 
a compact subset of the real line, and there are only a finite number of real eigenvalues. □ 

Remark 2. Note that, in the calculation of Haragu§ and Kapitula[21J the real eigenvalues play a slightly 
different role than other eigenvalues off of the imaginary axis. The fact that there are only a finite number 
of these indicates that there are only a finite number of values of the Floquet parameter for which there 
are real eigenvalues: n r (7) = for all but a finite number of values of the Floquet parameter 7 in their 
notation. 



2.1 Analysis of spec(d x £[u}) f| R 

We now move on to study the structure of spec(d x £[u]) (1 R more carefully. Suppose that fi e R. Clearly 
1 being an eigenvalue of M(/i) is a sufficient condition for fi G sj>ec(d x £[u\) and thus vanishing of 
D([i, 1) = a(fi) — a(—/j.) is a sufficient condition \x e spec(d x £[u]) n R. Notice that by the translation 
invariance of (7) we have D(0, 1) = by Nocthcr's theorem. The question is whether D(fi, 1) has any 
other real roots. If it does, then the eigenvalue problem (7) is spectrally unstable, due to the presence of 
a real non-zero element of spec(d x £[u}). In order to detect this instability, we calculate the orientation 
index 

I>(oo,l)D WM (0,l). 

As we will show, the negativity of this index is sufficient to imply a non- trivial intersection of spec(d x £[u]) 
with the real line. 

Lemma 1. The function D( ■ , 1) : R — > R is an odd function which satisfies the asymptotic relation 

lim D(u, 1) = =Foo. 

Proof. Clearly, D( ■ , 1) is an odd function of its argument, and hence it is sufficient to consider the limit 
as fi — > 00. To begin, define a new variable p = p}^T. Then from the asymptotic relations (10) we have 

a{p) = e-" + e -(-HV3iW2 + e -(-i-V3i)p/2 

a(p) = e P + e -(l+^)p/2 +e -(l-y3*)p/2 

where a(p) is the trace when you take p — ► —p. It follows that D(p, 1) = a(p) — a(p) behaves like — e p 
for large positive p, i.e. I «>0. This completes the proof. □ 
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From these results, we have the following theorem relating the sign of tr(M^ w (0)) to the stability 
of the underlying periodic wave. 

Theorem 1. Ifa"'(0) = tr( M MW1 (0)) > 0, then the number of roots of D(n, 1) (i.e. the number 
of periodic eigenvalues) on the positive real axis is odd. In particular spec(d x £[u}) nR* ^ I and the 
eigenvalue problem (7) is spectrally unstable. 

Proof. We show in Lemma 2 that D(0, 1) = D^O, 1) = D^(0, 1) = and D^(0, 1) = 2a'" (0). Thus, if 
a"'(0) > 0, then D(fi; 1) is positive for small positive values of \i. Since D(fi, 1) is negative for sufficiently 
large \i we know that D(±[i*,0) = for some /j,* 6 M \ {0}, which completes the proof. In the next 
section we establish the following formula for -0^^,(0, 1), the first non- vanishing derivative 



3 
E 



where again K is the classical action of the traveling wave ODE. Hence this "orientation index" can 
be expressed in terms of the Jacobian of the map between the constants of integration of the traveling 
wave ordinary differential equation (a,E;c) and the conserved quantities of the g-KdV (M,P,H). This 
orientation index is analogous to the quantity which is calculated in the stability theory of the solitary 
waves. □ 

It is important to notice the instability detected by Theorem 1 is an instability with respect to finite 
(bounded) wavelength perturbations. In the next section we will derive a modulational stability index 
which detects instability with respect to arbitrarily long wavelength perturbations. See the comments at 
the end of the section 3. The solitary wave solutions go unstable in the manner detected by Theorem 
1, through the creation of a pair of eigenvalues on the real axis. In general the periodic waves seem to 
first go unstable through the creation of a curve of spectrum which does not intersect the real axis, and 
later there is a secondary bifurcation resulting in a real eigenvalue. This phenomenon appears to have 
first been observed by Kapitula and Haragus,, who established that small amplitude periodic waves first 
go unstable at p = 2, as compared with p = 4 for the solitary waves. While we don't have a general 
proof of this we do show that, in the case of power law nonlinearity, there is a real periodic eigenvalue as 
well as a band of unstable spectrum connected to the origin. It is also worth noting that the analogous 
calculation for D(p,—1) shows that the number of anti-periodic eigenvalues on the real axis is always 
even. While this is not useful for proving the existence of instabilities it does eliminate some possible 
modes of instability. 



3 Local Analysis of the Period Map 

In this section, we turn our attention to studying the monodromy map M(/x) near the origin. To this 
end, we determine the asymptotic behavior of D(/i,e 1K ) as [i — > 0. We begin by proving that D(Q, e tK ) 
has a zero of multiplicity three at k = 0. It follows by directly computing the Jordan normal form of 
M(0) that A = 1 is an eigenvalue of algebraic multiplicity three and geometric multiplicity (generically) 
two. This fact reflects the following structure in the manifold of solutions to the ordinary differential 
equation defining the traveling waves: the traveling waves form a three parameter manifold, with traveling 
waves of constant period forming a two parameter submanifold. The two eigenvectors of the period map 
correspond to elements of the tangent plane to the submanifold of constant period solutions, while the 
third vector in the Jordan chain is associated to the normal to the constant period submanifold. 

Using perturbation theory appropriate to a Jordan block, as well as the Hamiltonian symmetry inher- 
ent in (7), we prove the three roots of D(fi, e tK ) bifurcate from /i = analytically in n in a neighborhood 
of K = 0, and derive a necessary and sufficient condition for modulational instability of the underlying 
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periodic wave u(x; a, E, c) in terms of derivatives of the monodromy operator. Note that this conclusion 
is somewhat unexpected: normally the eigenvalues of a non-trivial Jordan block do not bifurcate analyt- 
ically but instead admit a Puiseaux series in fractional powers. However because of the symmetries of 
the problem the admissible perturbations are severely restricted, resulting in a non-generic bifurcation. 

3.1 Calculation of the Period Map 

The first major calculation we present is an explicit calculation of the monodromy matrix at the origin 
in terms of the derivatives of the underlying periodic solution u with respect to the parameters. We do 
this by first computing a matrix valued solution to the ordinary differential equation satisfying the wrong 
initial condition: U(0, 0) is non-singular but not the identity. One can then multiply on the right by 
U _1 (0, 0) to find the monodromy matrix. We find that (as expected) the monodromy operator M(/i) 
has a non-trivial Jordan form when fi = 0. Our goal is then to utilize perturbation theory of Jordan 
blocks to calculate the normal form of the characteristic polynomial in a neighborhood of /i = 0, A = 1, 
where A is the eigenvalue parameter of the monodromy operator. 

To begin we write the above third order eigenvalue problem as a first order system as in (8). In 
particular, notice that tr( H(x; /i)) = for all and thus det($(x;/x)) = 1 for all /ieC, implying 
det(M(^)) = 1. In order to calculate a matrix solution <!>(x;^), we must first find three linearly inde- 
pendent solutions of the above system. In general, this is a daunting task, but since the above system 
with /I = arises as the Frechet derivative (linearization) of an integrable ordinary differential equation 
this can be done by considering infinitesimal variations of the constants of integration in the defining 
ordinary differential equation, and thus generating the tangent space. As noted earlier the solutions 
u(x — x ; a, E, c) constitute a 4-dimensional solution manifold of (1) parameterized by xo, a, E, c. The so- 
lutions of the linearized operator space are given by the generators -4- , 4- , and -4^ acting on the solution 
u(x; a, E, c). The action of the generator 4- c is somewhat different and is connected with the generalized 
null-space. This will become important in the next section. 

Proposition 2. Let u(x; a, E, c) be the solution to the traveling wave equation (5) satisfying u(0; a, E, c) = 
u^,u x (0; a, E, c) = 0. A basis of solutions to the third order system 

Y x = H(x; 0)Y 

is given by 

{u X: U XX: u xxx ) 
(u a: U aX7 U aX x) 

{u E ,u Ex ,u Exx ). 
A particular solution to the inhomogeneous problem 

Y x = H(x; 0)Y + W 

where W l — (0,0, u x ) is given by 

Y^ = (u Cl u CXl u cxx ). 

Proof. A straightforward calculation. Notice that it follows that d x C[u] (— u c ) = u x . □ 

The fact that u a , u E are not periodic - they exhibit secular growth due to the variation of the period 
with respect to the parameters - gives an indication that the eigenspaces of the monodromy at \i — are 
not semi-simple, and hence we expect the existence of a non-trivial Jordan block of the monodromy map 
M(0). 



Yl = 

n = 
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By the above proposition, three linearly independent solutions of (8) corresponding to fi = are given 

by 

(u'(x;a,E,c) \ / u a (x;a,E,c) \ / u E (x;a, E,c) \ 

u"(x;a,E,c) Y 2 (x) = u' a (x;a,E,c) Y 3 (x) = U ' s (x;a,£,c) . (11) 
u'"(a;;a,£,c) / \ u'^(x; a, E, c) / \ u" E (x;a,E,c) ) 

By hypothesis, for any a, £7, c € R the solution u satisfies 

u(0;a,E,c) = u_ = u{T;a,E,c) (12) 
d x u(0; a, £, c) = = d x u(T;a,E,c) (13) 
d xx u(0;a,E, c) = a-/(u_) + cu_ = d xx u(T; a, E, c). (14) 

Moreover, from equation (1) it follows that 

■u raa; (0; a, £, c) = cu x (0; a, E,c)--^ {f(u(x; a, E, c))) | x=0 = 0. 

Defining U(x, 0) = ^2(2^), ^3(2;)] to be the corresponding solution matrix, then direct calculations 

yield 

(0 d a u- d E u- \ 

a- f(u-) + cu- (15) 

\ + {c- f'{u-))d a u- (c- f'(u-))d E u- J 

Note that differentiating the relation _E — V(m_) = gives the relation — V {u^)8eU- = dct(U(0, 0)) = 
— 1, so these solutions are linearly independent at x = 0, and hence for all x. Thus we can compute 
U(T, 0) and right-multiply by U~ 1 (0, 0) to give the monodromy M(0) 

The matrix U(T, 0) can be calculated by differentiating (12)-(14) with respect to the parameters 
a and E by use of the chain rule. For example, differentiating the relation (12) with respect to the 
parameter E gives 

d E u{T; a, c, E) + — (T; a, E, c) T E (a, c, E) = 8 E u_ . 

Since the derivative vanishes at the period points this implies d E u(T) = d E U-. Continuing in this manner 
gives the following expression for the change in tis matrix solution across the period: 



U(T,0) = U(0,0)+ I y'(w_;a,c)T a V'(u-;a,c)T E |. (](,) 


In particular, we find that U(T, 0) — U(0, 0) is a rank one matrix, which naturally leads to the following 
proposition. 

Proposition 3. There exists a basis in R 3 such that the monodromy matrix M(p) evaluated at fi = 
takes the following Jordan normal form: 




M(0) ~ 1 a (17) 



where a 7^ as long as T a and T E do not simultaneously vanish. In particular, the monodromy operator 
at = has a single eigenvalue of A = 1 with algebraic multiplicity three and geometric multiplicity two 
as long as the period is not at a critical point with respect to the parameters a, E for fixed wavespeed c. 
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Proof. Recall det( U(0, 0)) = -V'{u-)d E u- = -1, so U(0,0) is invertible. Multiplying the above 
expression on the right by the matrix U~ 1 (0, 0) yields the monodromy matrix at the origin 

MO = 0) := / + w® v U _1 (0) 

where w = (0, 1,0) T and v = (0, V'(uJ)T a , V'{u_)T e ) t . Next, notice that 

/ -T a -T E \ 
U(0,0) = w®v 
V / 

and hence defining N := U _1 (0) M(0) U(0) gives the equation 

/ 1 -T a -T E 
N= 1 
\ 1 

It follows that 

Ker{ N - I) = span {{I, 0, 0) T , (0, T E , -T a f} 

Now, take Vs := (0, — T a , —T E ) and notice that v% £ Ker(N — I). The Jordan structure then follows by 
noticing then that ( N - I)v 3 = (T 2 + T§)(1, 0, 0) T e ker(N - I). 

For traveling waves below the separatrix a result of Schaaf (see appendix) shows that 2 the period is 
a strictly increasing function of the energy, T E > and thus the gencricity condition is always met. In 
other situations we will assume that this condition is met unless otherwise stated. □ 



3.2 Asymptotic Analysis of D(fj,, n) near (/i, k) = (0, 0) 

We now analyze the characteristic polynomial of M(/i) in a neighborhood of /x = by considering M(/i) 
as a small perturbation of the matrix M(0) constructed above. It is well understood how the eigenvalues 
of a Jordan block bifurcate under perturbation: see Kato[25] or Moro, Burke and Overton[30]. It is worth 
noting, however, that in this case the bifurcation is highly non-generic due to the constraints imposed by 
the symmetry of the problem. 

Recall from Proposition 1 that the spectrum near /i = is continuous. By the analyticity of M(/i) 
in a neighborhood of fi = 0, we can expand M(/x) for small \x as 

M(u) = M(0) + uM^O) + fL. M w (0) + O (ImI 3 ) 

where M^O) = [AffV] and M w (0) = [M$]. If one makes a similarity transform M(/i) = V" 1 M(it)V 
so that M(0) is in the Jordan normal form (17) then a direct calculation using the above second order 
expansion of M(/i) implies that in a neighborhood of fj, = 0, the characteristic polynomial can be 
expressed as 

D(u, e lK ) = det ((M(ii) - 7) - {e lK - l)l) 

= -i) 3 + J | 2 ^tr(M (1 (0))+^-tr(M fJfJ (0))j 
- V ^M^a + f, 2 Q (tr(M,)) 2 - ^tr(Mj) - |m 3 (2 2 } ^ 

+ fi 3 (det (M M (0)) +<tS) +0(4), (18) 

2 Under some mild assumptions on the nonlinearity, which are satisfied for the power law nonlinearity. 
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where rj — e lK — 1, S represents mixed terms from M„(0) and M w (0), a is as in Proposition 3, and the 
notation O (4) represents terms whose degree is four or higher. Notice there are no other /j, 3 terms since 
M(0) — / has rank one. Our next goal is to determine the dominant balance of the equation D(fi, e %K ) = 
in a neighborhood of (/i, n) = (0, 0). 

A useful construction for implicit function calculations of this type is that of the Newton diagram, 
which is a subset of the non-negative integer lattice. A vertex is included if the coefficient of rf^ 1 ^ 
in (18) is non-zero, otherwise the vertex is not included. The lower convex hull of the Newton diagram 
is made up of a collection of line segments. For each line segment of the lower convex hull let m be the 
horizontal length of the segment and s the slope of the segment. Corresponding to each such line segment 
there are m distinct solution branches of the form 

^) = E a *V, 

i 

where a[ k ^ ^ and k ranges from 1 to m. For details see the book of Baumgartel[3] or Hilton[22]. This 
is equivalent to the method of "dominant balance" presented in textbooks on asymptotic methods but 
is somewhat more systematic. For instance, in our case if the coefficient of the rj 1 ^ 1 term (-aM^) is 
non- vanishing then there are two solution branches in which r\ has an expansion in powers of fii and one 
with an expansion in integer powers. These correspond to the breaking up of the 2x2 and lxl Jordan 
blocks respectively. However as mentioned above the symmetry M(/i) <~ M _1 (— /i) causes a number of 
terms in (18) to vanish, which leads to an expansion in integer powers of \x. This is the content of the 
next lemma. 

Lemma 2. The equation D(/j,, e lK ) = has the following normal form in a neighborhood of (fi, k) = (0,0): 

_ (i «)3 + ^ tr ( M w (0)) + ^ tr ( M^(0)) + 0(4) = 
whose Newton diagram is depicted in Figure 2. 

Proof. Define functions a, b, c on a neighborhood of fi — by 

det[( M(/x) -I) - (e lK - 1)1] = -if + (o(/i) - 3)r/ 2 + 6(^)77 + c(p). (19) 

where r\ = e %K — 1. Notice in particular that r\ — in + O (k 2 ) in a neighborhood of k = 0. By (18), it 
follows 

2 3 

a(p) = tr(M( A1 )) = 3 + / itr(M AJ (0)) + ^-tr(M w (0)) + ^tr(M MW (0)) + O(| A1 | 4 ) 

KM) = \ ( tr(( M(jm) -If)- tr( Mfjj) I) 2 ) = -/xAfga - M 2 Q tr( M,) 2 - \ tr( M 2 ) - |m 3 , 2 ) + 0{\tf) 

cQj) = det(M(/i) - I) = -a(MWM$ - M^M^)n 2 + (det(M /1 (0)) + aS) + G(| M | 4 ) 

Using the symmetry M(— fi) ~ M(^i) -1 , we have 

c(-/x) = det[M(-/x)-i] 

= detain)]- 1 det[I - M(n)] 
= -det[M(/x)-i] 

= - c (m) 

since det[M(/i)] = 1 for all \i e C. Hence c is an odd function of [i. Also, since M(0) — I has rank 
one, (18) along with the above analysis implies that c(p) = C(|/i| 3 ), from which it follows c"(0) = 
-2*(M^M$ - M$M[% = 0. 
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Figure 1: The Newton diagram corresponding to the asymptotic expansion of D{fi, e lK ) = in a neigh- 
borhood of (/U, k) = (0, 0) is shown to 0(|/x| 3 ). Terms associated to open circles are shown to vanish due 
to the natural symmetries inherent in (1). The grey circles are non- vanishing terms which are a part of 
the lower convex hull. The black circles lie above the lower convex hull and thus do not contribute to the 
leading order asymptotics. 



Similarly, using (19) we have 



det[M(/i) - AI] 



A 3 det M(-p) 




+ a(-fj.) 



= -(A - l) 3 - a(- M )A(A - l) 2 + 6(- M )A 2 (A - 1) - 

Comparing the A 2 and A 3 terms above with those in (19) we get the relations 

b(fi) — 2a(/i) — a(— fj) — 3, and 
a(p) - b(p) + c(/x) = 3. 




Since c(/z) = C(|m| 3 ), these relations imply tr(M M (0)) = a'(0) = and -crM$ = b'(0) = 0. By 
recalling a ^ from 3, this implies M§\ — 0. Moreover, we know that o"(0) = b"(0) and hence 
6"(0) = tr(M w (0)). Also, we have b"'(6) = 3a'" (0) and c"'(0) = b"'(0) - a"'(0) = 2a"' (0) and hence 
c"'(0) = 2tr(M jUjUjU (0)). The corollary follows by analyzing equation (18) as well as the corresponding 
Newton diagram (see Figure 1). □ 

From this it follows that, in the neighborhood of the origin, the leading order piece of the periodic 
Evans function is a homogeneous cubic polynomial in K,fi. The implicit function theorem fails, but in a 
trivial way that is easily corrected, leading to the following theorem: 



Theorem 2. With the above notation, define 



A(/;«) = i(tr(M w (0))) 3 - 3 ( tr( M^O))) 2 , 



(20) 



where f denotes the dependence on the non-linearity used in (1), and suppose tr(M^ w (0)) 7^ 0. If A > 0, 
then the imaginary axis in the neighborhood in the origin is in the spectrum with multiplicity three. If 
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Figure 2: When A(/; u) < 0, the local normal form of spec(d x £[u]) consists of a segment of the imaginary 
axis union with two straight lines making equal angles with the imaginary axis. Notice that these lines 
intersect at the origin, corresponding to the fact that 1 is an eigenvalue of M(0) with algebraic multiplicity 
three. 

A < then the imaginary axis in a neighborhood of the origin is in the spectrum with multiplicity one, 
together with two curves which are tangent to lines through the origin with angle arg(ij/2,3) - see Figure 
2. In particular in the latter case the periodic wave is modulationally unstable. 

Proof. Since the leading order piece of the Evan's function is homogeneous it suggests working with a 
projective coordinate y = Making such a change of variables leads to the equation 

2 3 

1 - V - tr ( M„„(0)) + V - tr ( M WM (0)) + kE[k, y) = 

where E(k, y) is continuous in a neighborhood of the origin. The implicit function theorem applies in a 
neighborhood of (y = yi.2,3, « = 0) as long as the roots 2/1.2. 3 of the above cubic in y are distinct, which 
is true as long as the discriminant A is not zero. In terms of the original variable /x we have the three 
solution branches 

A»l,2,3 = iyi,2,3K + 0(k 2 ) 

This cubic has three real roots when A > 0, giving three branches of spectrum emerging from the origin 
tangent to the imaginary axis. It is clear from symmetry that these must in fact lie on the imaginary 
axis, giving a interval of spectrum of multiplicity three along the imaginary axis. In the case that the 
discriminant is negative there is one real root and two complex conjugate roots, giving one branch of 
spectrum along the imaginary axis and two branches emerging in the complex plane. □ 

Remark 3. First we note that tr(M jUjU (0)) < is a sufficient condition for modulational instability of 
the periodic wave. 

Secondly we note that the Newton diagram is independent of tr(M MM (0)) but the dominant balance 
changes if tr(M ([tjUjU (0)) should happen to vanish. Later we will give a formula for tr(M ([tjUjU (0)) in terms 
of the Jacobian of a map, and we will see that vanishing of this quantity signals a change in the Jordan 
structure of the underlying linearized operator. 

Finally we note that this agrees with the result of Bottman and Deconinck[8], in which they consid- 
ered cnoidal wave solutions to the KdV. Using the algebro- geometric techniques of Belokolos, Bobenko, 



3 LOCAL ANALYSIS OF THE PERIOD MAP 



16 



Enol'skii, Its and Matveev[5] they explicitly computed the spectrum of the linearized operator and found 
that such solutions are always linearly stable. Their results prove that an interval of the imaginary axis 
containing the origin is in the spectrum of the linearized operator, with multiplicity three. Our results 
imply this is a generic phenomenon: either one has an interval of spectrum of multiplicity three about 
the origin, or one has three curves intersecting at the origin. In the KdV case the discriminant A is 
expressible in terms of elliptic functions in this case and must be positive, although we have not tried to 
show this. 

It is instructive to compare this theorem with that of Theorem 1. By Theorem 2, the sign of 
tr( M /JMM (0)) has no effect on the spectral stability of the underlying periodic wave in a sufficiently small 
neighborhood of the origin. However, Theorem 1 guarantees the existence of unstable real spectrum if 
sgn(tr(M MAIAI (0))) > 0. To reconcile these results, notice Proposition (1) implies there is no unstable 
real spectrum sufficiently close to the origin. Thus, the instability brought on by tr( M MW (0)) > is not 
local to /i = 0, and hence should not be detected by the quantity A(a, E, c). 

Our next goal is to try to evaluate the modulational stability index A(a, E, c) as well as the finite 
wavelength orientation index tr( M AtMM (0)) in terms of the conserved quantities of (1). This can be 
done very explicitly. Notice that while we have chosen to express the coefficients as tr( M A , AI (0)) and 
tr( M /JAIAI (0)), which suggests that they arise at second order and third order in a perturbation calculation 
for small /it, these quantities can be expressed in terms of quantities which arise at first and second order 
in fi due to the invariance of the problem under the map /j — > — /i. Further while all of the first order 
terms contribute only a few terms which are second order actually contribute - these are the terms which 
are associated to the minors of the off-diagonal piece of the unperturbed Jordan form. These second 
order terms are explicitly computable via a single quadrature. 

Theorem 3. We have the following identities: 



where T, M, P are the period, mass, and momentum of the underlying traveling wave and a, E, c parame- 
terize the family of traveling waves. Thus the modulational stability index has the following representation 



Proof. Let Wi(x;/T), i = 1,2,3, be three linearly independent solutions of (1), and let W(x,/i) be the 
solution matrix with columns Wi. Expanding the above solutions in powers of /i as 



Using Proposition 2, we choose Wi{x) = Yi(x) where the vectors Yi(x) are defined in equation (11). The 
higher order terms in the above expansion yield 



tr(M w )|„ =0 = {T,P}B,c + 2{M,P} o ,B 
tr(M WM )| M=0 = -^{T,M,P} U}E}C 




Wi{x,n) = w°{x) + »wj{x) + ^ 2 wj(x) + 0(\fi\ 3 ) 



and substituting them into (7), the leading order equation becomes 



-£«;?(*)= Hfo <))«;?(*). 



d_ 

dx 



wj{x) = H(x;0)wl(x) + Vt 1 



(x), j>l 




where Vf = [0,0, — (u^ )ij and (v)i denotes the first component of the vector v. Notice that for 
each of the higher order terms j > 1, we require w l j(0) = 0. This implies that W(0,/x) = U(0, 0) in 
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a neighborhood of /i = 0, where U(0, 0) is defined in (15). In the case j = 1, the i = 1 equation is 
equivalent to the equation L w\ = u x . It follows again from Proposition 2 that we can choose 



w\{x) = 




u E 



Notice the above coefficients of Y 2 (x) and Y 3 (x) are determined by differentiating E — V(u-; a, c) = with 
respect to the parameters a, E, and c. Moreover, using variation of parameters as well as the identities 
{u x ,u} x> e = -1 and {u,u x } x . a = u, we choose 



W(a;,0) / W(z,0)- 1 V^ 1 (z)dz 
Jo 




* Io( W i~ rf Z _ U a JiT ( W i ~ ) 1 dz + U E Jq (wj ~ 1 ) 1 U dz \ 

Io(wj~ 1 )i{u,u x } at E dz- u ax J*(wt~ 1 ) 1 dz + u Ex f*(w : l~ 1 ) 1 udz (21) 
fxxxfa(wj~ 1 )i{u,u x } a , E dz ~u axx ]*(w{~ 1 )idz + u Exx $ X {w{~ 1 ) l udz j 

for w\, i — 2,3, and for i > 1. Finally, by (16), we have the following expression valid as /i — > 0: 

/ 0(H 2 ) 2 0(H) 0(H) \ 

-5W( M )= M ^'(w-)(-T c + M _T a -4r £ )+0(|M| 2 ) V'(u-)T a + <D(\n\) V'{uJ)T E + 0(H) , 

V " 0(H 2 ) ' 0(H) 0(H) / 

where 5W(/i) := ~W(x, /u) |^_ , and the C(/i) and C(/i 2 ) terms above are computed using (21). Note 
that all of the (D([i,fi 2 ) terms in the above are necessary for the calculation, however we do not write 
them out. Recalling that our choice of basis implies det(W(0, fj,)) = —1, we have 

D(p,l) = -det(5W(/x)) 

= -i{T,M,P} a ^ cM 3 + G(| M | 4 ), 

from which the expression for tr(M Wft (0)) follows by Theorem 1. Moreover, it follows from Lemma 2 
and the fact M(/i) = iSW(ju) W(0, 0) _1 + I and a rather tedious calculation that 

tr(M w (0)) = -2 / i- 2 tr(cof(M( M )-/))| A(=0 
= -2^{T,P} E . c -{M,P} E . a 

as claimed. □ 

Corollary 1. {T, M, P} a ,E,c < is a sufficient condition for a non-trivial intersection of spec(d x £[u]) 
with the real axis. 

Proof. This is now clear from Theorems 1 and 3. □ 

At this point we can make a connection to the stability theory for the solitary waves. 

Corollary 2. In the case of power-law nonlinearity and wavespeed c > 0, there are always unstable 
periodic traveling waves in a neighborhood of the solitary wave (a = E = 0) if p > 4. Moreover, such long 
wavelength periodic waves exhibit a modulational instability if and only if p > 4. 



3 LOCAL ANALYSIS OF THE PERIOD MAP 



18 



Proof. First, note that the scaling invariance in equation (4) implies the periodic solution u(x; a, E, c) 
satisfies 



u(x; a, E, c) = c^u (c^x; lj 



which allows us to compute T c , M c , P c explicitly as follows: 

-J-n ±E 







a 


T c = 


-It 






2c 










M c = 








\pc 










Pc = 


i 1 - 




\pc 





pc pc 



■ 1 r» / v w w / 

2c J pc pc 

T\ 2 a(p+l) , 2 v E(2+p) fn 2 , 

T c + — ]u_ (P a - T a u_) (P E - T E u_) , 

2c J pc x ' pc x ' 

where T c follows from equation (6). Since we know that Pe = 2T C and P a — 2M C the above serves 
to simplify the last row and column. When a and E are small there are two turning points ri,r 2 in 
the neighborhood of the origin and a third turning point r 3 which is bounded away from the origin. In 
the solitary wave limit o,J? we have n — r 2 = 0(Va 2 — 2E). In this limit we have the following 
asymptotics for small (a, E) 

M(a,E,l) = 0(1) 
P(a,E,l) = 0(l) 
T(a, E, 1) = O (ln(a 2 - 2E)) 

T a (a,E 1 l) = o(- T ^—) =M E 



T E {a,E 1 l) = 



o?-2E 
1 



S a 2 -2E / 

Thus the asymptotically largest minor of {T 7 M,P} a E c is —T E M a P c , from which it follows 

{T, M, P} a , E ,c - -T E M a (---p\P 



pc 2c 

as a, E — > 0. It is known that for traveling waves below the separatrix (see appendix) that under 
some minor convexity assumptions Te > 0. It can also be shown that (see appendix) M o (a,0) < for 
E = and a sufficiently small. Thus the orientation index {T, M, P} a ,E,c is negative for p > 4 and a, E 
sufficiently small (in other words sufficiently close to the solitary wave) and positive for p < 4 and a, E 
sufficiently small. This also follows, of course, from Gardner's long-wavelength theory[17] but it provides 
a good check for the present theory. 

To prove the second claim, notice the above asymptotics implies 

tr(M w( 0)) ~ T £ (l-1 

in the limit as a, E — > 0. Hence, it follows that sgn A(a, E, c) = sgn(4 — p) for a, E sufficiently small, and 
periodic waves of sufficiently long period are also modulationally unstable for p > 4. 

□ 

Remark 4. It is worth noting that the instability mechanism detected by the discriminant A is not 
present in the solitary wave case: in the solitary wave limit the bands of spectrum connected to the origin 
collapse to the origin. This instability does not appear to follow from Gardner's calculation: Gardner 
shows that the point eigenvalue of the solitary wave opens into a small loop of spectrum, predicting the 
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real eigenvalues detected by {T, M, P} a ^E,c, but the modulational instability detected by A is not detected. 
Thus suggests the heuristic that periodic solutions should go unstable before the solitary waves. The small 
amplitude stability calculation of Haragus and Kapitula for the generalized KdV equation amounts to a 
calculation of this discriminant in that limiting case, and their proof that the small amplitude waves go 
unstable at p = 2 is the first result we are aware of along these lines. 

We believe that a small amplitude analysis of A(a, E, c) should be possible. It follows by a simple 
calculation that A = at the stationary solution. By expanding near by solutions in terms of amplitude 
instead of the energy E, we believe the first non-zero term of the discriminant should be proportional to 
a polynomial which switches signs at p = 2, thus recovering the small amplitude result of Haragu§ and 
Kapitula [21] . We have not as yet carried out such an analysis. 

Using the identities derived in Appendix 1, we now have a sufficient criterion for the existence of a 
non-trivial intersection of spec(d x £[u\) with the real axis in terms of the conserved quantities M , P and 
H of the gKdV flow, as well as a necessary and sufficient condition for understanding the normal form of 
the spectrum in a neighborhood of the origin. It is a rather striking fact that both of these indices can 
be expressed entirely in terms of the conserved quantities of the flow. The monodromy itself depends on 
U-(a, E, c), the classical turning point of the traveling wave, as well as various functions and derivatives 
of this quantity, but the indices themselves only depend on the conserved quantities. This is, of course, 
the Whitham philosophy, but we are only aware of a few cases (other than the integrable calculations, 
which are very special) in which make this rigorous. 

In the next section we outline the connections of this calculation to a calculation based more directly 
on the linearized operator. While not strictly necessary this calculation is useful since it clarifies the 
way in which various bifurcations can occur. In this section we calculate the null-space and generalized 
null-space of the linearized operator and sketch a perturbation calculation analogous to the one given for 
the Evan's function. 

4 Local Analysis of spec(d x £[u]) via the Floquet-Boch Decompo- 
sition 

4.1 Floquet-Bloch Decomposition 

In this section we sketch an approach to this problem working directly with the linearized operator rather 
than with the Evan's function. While these two approaches are presumably equivalent the former seems 
less straightforward than the latter. In particular it is not clear how one might derive the orientation 
index in this way, and the calculation of the modulational stability index gives a quantity which seems 
much less transparent. Nevertheless we present an outline of this calculation (omitting some details) 
since it does give some insight into the results of the previous section. 

From Floquet theory, we know any bounded eigenfunction v(x) of c^Z^u] must satisfy 

v(x +T) = e^v{x) 

for some 7 e [— 7r,7r]. The quantity e n is known as the Floquet multiplier of the eigenfunction v. It 
follows any eigenfunction v(x) can be represented in the form v(x) = e %lx l T P[x) where P(x + T) = P{x). 
The fact that d x C[u]v(x) = fiv(x) for some \i <G C implies 

e l7x/T 'J 7 r 7 [u]P(x) = ^ x l T P{x) 

where J 7 = (d x + i^) and £ 7 [u] = — (d x + i^) 2 — f'{u) + c arc the so-called Bloch operators. This 
suggests fixing a 7 <G [— n, tt] and considering the eigenvalue problem for the operator J 7 £ 7 [w] on the 
Hilbert space H = L 2 (R/TZ;C). This procedure is known as a Bloch decomposition of the eigenvalue 
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problem (7) and we consider the Bloch operators operators as acting on Ti. Notice for 7 / the 
operators J 7 £ 7 [w] are closed in this space with compactly embedded domain H 3 (R/TZ;C). It follows 
that these operators have compact resolvent and hence their spectra consists of only point spectra with 
finite algebraic multiplicities. Moreover, one has 

spec(d x £[u]) = [J spec (J 7 £ 7 [u]) . 

7£[— 7T,7r] 

Thus, this decomposition reduces the problem of locating the continuous spectrum of the operator d x C[u] 
on L 2 to the problem of determining the discrete spectrum of a one parameter family of operators 
{ J 7 £ 7 [tt]} 7e [_ 7r . x ] on H. Our first goal is to understand the nature of the spectrum of the operator 
J £oM at the origin. Notice in particular that for 71,72 € [— n, n], J 71 £ 71 [u] is a compact perturbation 
of J 72 £ 72 [it], and hence routine calculations prove the above parameterization of spec(J 7 £ 7 [u]) is in fact 
continuous. We then consider the operator J 1 L 1 [u] for I7I <C 1, treating it as a small perturbation of 
Jq£o M ! w hh our end goal being to study how the spectrum bifurcates from the 7 = state. 

We begin with analyzing the generalized periodic null space of the operator d x C[u], denoted N g {d x C[u\) - 

uZi N((d x c[u]r). 

Proposition 4. Suppose that the Jacobians {T, M} a _ E and {T, P} a ,E do not simultaneously vanish, and 
{T, M, P} a ,E,c 7^ 0. Then zero is an eigenvalue of the operator d x C[u] = J £ [u] considered on H of 
algebraic multiplicity three and geometric multiplicity two. In the case {T, M} a _E ^ 0, the functions 

4> = {T 7 u} at E, ^0 = 1, 

4>! = {T, M} at E U x , tpl= 4>2(s)ds, 

Jo 

2 = { Ul T,M} a . E , c V 2 = {T,M} B , c + {T,M} a , EU , 

provide a basis for N g {d x C[u}) and N g (C[u]d x ) respectively. Specifically we have the relations 

d x £[u](fio = £[u]d x ip = 

d x C[u]4>i = £[u]d x ibi = -ip 2 

d x £[u\4> 2 = -4>i C[u]d x %l> 2 = 0. 

Proof. The constants above are chosen for convenience, and the functions above are not normalized. For 
instance, 4>2 can be any multiple of u x and similarly ipo any constant. Also, the ordering is chosen so 
that (<f)j,ipk) = for i 7^ k. Notice this proposition does not follow directly from Proposition 2 since the 
functions u ai ue and u c are not in general T-periodic, and one must chose linear combinations which are 
periodic and thus belong to Ti. 

First observe that (16) implies (j)o and <pi are T-periodic and belong to iV(J 7 £ 7 [u]). In particular, 
Corollary 3 from the appendix implies N(d x C[u\) = spanj^o, 4>i} when considered as an operator on Ti. 
The fact that the monodromy at the origin is the identity plus a rank one perturbation suggests that 
there are two linear combinations that can be chosen to be periodic. Specifically we define 



u a T a u a dx 
ue T e Jo u E dx 
u c T c J Q u c dx 



= {u,T,M} a , E ,c 



and it is clear from (16) that <f)2 € H and Jo£oM02 = 4>i as claimed. Thus, if {T,M} at E 7^ 0, <f>2 gives a 
function in N((J C [u}) 2 ) - N(J C [u}). 

Similarly, Corollary 3 implies that V'o and ^2 arc belong to N(£[u]d x ), and are linearly independent 
provided that {T, M} a ^E 7^ 0. Moreover, a it is clear from construction that ipi € H and a straightforward 
computation shows that ib\ belongs to N((£ [u]J ) 2 ) — N(£ [ u ]Jo) as claimed. 



4 LOCAL ANALYSIS OF SPEC(<9 x £[t/]) VIA THE FLOQUET-BOCH DECOMPOSITION 



21 



In order to complete the proof, we must now show these three functions comprise the entire generalized 
null space of J £o[ u ] on To this end, we prove that neither of the functions </> , 4>2 belong to the range 
of J £o[u] by appealing to the Fredholm alternative. It follows that the equation J C [u]v — </> has a 
solution in Ji if and only if the following solvability conditions are simultaneously satisfied: 

(l,0o) = {T,M} o>B = 0, and 

(U^O) = \{T,P} a ,E = 0. 

Thus, if either {T,M} aE or {T, P} a ,E are non-zero, then N ((J Q C Q [u]) 2 ) — N(J £ [u]) = span{0 2 }- 
Similarly, N ((J £o[ u ]) 3 ) — N ((J n £ [ u }) 2 ) 7^ if and only if the equation L n v = (f> 2 has a solution in H, 
i.e. if and only if 

(u.fo) = l{T,M,P} a , B , c = 0, 

which finishes the proof. 

A similar construction in the case {T, M} a ^ E = but {T, P} a ,E 7^ gives a basis in this case. □ 

Remark 5. It is worth remarking in some detail on the physical significance of these conditions and the 
relationship to the Whitham modulation theory. Obviously (a, E, c) are constants of integration arising 
in the ordinary differential equation defining the traveling wave, and T, M, P are constants of the PDE 
evolution. One of the main ideas of the Whitham theory is to locally parameterize the wave by the 
constants of motion. The non-vanishing of the Jacobians is exactly what allows one to do this. Non- 
vanishing of {T, M, P} a ,E,c is equivalent to demanding that locally the map (a, E, c) 1— » (T, M, P) have 
a unique C 1 inverse - in other words the conserved quantities (T,M,P) are good local coordinates for 
the family of traveling waves. Similarly non-vanishing of one of {T, M} a ,E and {T, P} a ,E is, at least for 
periodic waves below the separatrix, equivalent to demanding that the matrix 

f T a M a P a \ 
\ T E M e Pe ) 

have full rank, which is equivalent to demanding that the map (a,E) (T,M,P) (at fixed c) have a 
unique C 1 inverse - in other words two of the conserved quantities give a smooth parameterization of 
the family of traveling waves of fixed wave-speed. As long as E ^ we can use the identities developed 
in the appendix to eliminate T in favor of H. Thus in the case £ ^ (which does not include the 
solitary wave wave) the null-space being two dimensional is equivalent to two of the conserved quantities 
(M, P, H ) giving a C 1 parameterization of the traveling wave solutions at constant wavespeed, and the 
space N((JqCo[u}) 2 ) — N(JoCq[u]) being one dimensional is equivalent to the three conserved quantities 
(M, P, H) giving a C 1 parameterization of the full family of traveling waves. 

Notice it follows the vanishing of {T, M, P} a ,E,c is connected with a change in the Jordan structure 
of the linearized operator Jo£q[u] considered on TL: {T, M, P} a ,E,c 7^ ensures the existence of a non- 
trivial Jordan piece in the generalized null space of of dimension exactly one. Moreover, it guarantees 
that the variations in the constants associated to the family of traveling wave solutions by reducing (1) 
to quadrature are enough to generate the entire generalized periodic null space of the operator JqCq [u]. 
Henceforth, we shall assume {T, M, P} a ,E,c 7^ and that {T, M} a ,E 7^ - trivial modifications are 
necessary if {T,M} a ,E vanishes but {T, P} a ,E does not. 

4.2 Analyticity of Eigenvalues Bifurcating from n — 

Our next goal then is to consider the operator J 7 £ 7 [u] for small 7, treating it as a small perturbation 
of J £ [u]. To this end, notice that if we define L := Jq£q{u], L\ := £o[ u ] — 2<9^, and L 2 := — 3d x , it 
follows that 

J 7 £ 7 [u] = L + eL\ + e 2 L 2 — e 3 , 
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where e is related to the Floquet exponent via s = By Proposition 4, we know the operator Lq has 
three periodic eigenvalues at the origin. Our present goal is to determine how these eigenvalues bifurcate 
from the 7 = state. In this section we only sketch the relevant details - for similar calculations see the 
papers of Ivey and Lafortune[23], or Kapitula, Kutz and Sanstcdc.[24] 

Since the Hilbert space TL consists of T-periodic functions, eigenvalues of J 7 £ 7 [w] correspond to 1 
being an eigenvalue of the monodromy operator <&(T; /j,, e) for to the eigenvalue problem J 7 £ 7 [u]u = 
Thus, it is natural to introduce the following "modified" periodic Evans function 

D (ji,e) = det(&{T;ii,e)-I). 

Notice that D (fj,, e) is clearly an analytic function of the two complex variables fj, and e. Our first goal 
then is to analyze the possible behavior of of the solutions of D (n,e) = in a small neighborhood of 
(0,0). 

Lemma 3. Let F(x, y) be a complex valued function of two complex variables x and y which is analytic 
in a neighborhood of (0, 0) E C 2 . Moreover, suppose that F(0, 0) = ^(0, 0) = F xx (0, 0) = 0, F XXX (Q) 7^ 0, 
and ^(0,0) = 0. Then for small y, the equation F(x,y) = has three roots in a neighborhood of the 
origin. Moreover, these roots are given by (x,y) = (fj(y),y), j = 1,2,3, where the fj satisfy one of the 
following conditions: 

(i) One function fj can be expressed as a Puiseux series as fj(y) — Y^?=i a ny n ^ 2 ^ n a neighborhood of 
y = 0, where a\ 7^ 0. 

(ii) Two of the functions fj admits a Puiseux series representation of the form fj(y) — Yl^=2 a nV n ^ 
in a neighborhood ofy — 0, where a 2 7^ 0. 

(Hi) All three functions fj are O(e) and are analytic in y in a neighborhood of y — 0, i.e. they can be 
represented as fj(y) = J2™=i a nU n where a\ 7^ assuming F yyy (0, 0) 7^ 0. 

In the case (Hi), if F yyy (0,0) = then all three eigenvalues are analytic in e, with two eigenvalues of 
order 0(\e\) and the remaining eigenvalue of order at least 0(\e\ 2 ). 

Proof. By the Weierstrass preparation theorem, the function F(x,y) can be expressed as 

F(x,y) = (x 3 + T] 2 (y)x 2 + n 1 (y)x + r] (y)) h(x,y) 

for small x and y, where each r]j is analytic, and h is analytic satisfying h(0, 0) 7^ 0. It follows the three 
roots of F(x,y) near (0,0) are determine by the cubic polynomial G(x, y) = x 3 + n 2 (y)x 2 +rji{y)x + rj {y). 
By hypothesis, we have that f]j(0) = for j = 0, 1,2, T7q(0) = 0, and t]q'(0) 7^ 0. Hence, the Newton 
diagram for the equation G(x, y) = is the same as that in figure 1, from which the lemma follows. □ 

We now wish to apply Lemma 3 to the equation Dq(/i, e) = 0, with x = \x and y = e, and use the Fred- 
holm alternative to show only possibility (Hi) can occur. Notice that Theorem 3 implies -^-^ Dq(/j,; 0) = 

for k = 0,1,2 and, moreover, -^Dq{^l-.Q) 7^ under the assumption {T, M, P} a ,E,c 7^ 0. To apply 
Lemma 3 then, we need the following lemma. 

Lemma 4. We have ^D (0,0) = 0. 

Proof. This proof proceeds much like that of Theorem 3. Defining W(x;/j,,e) to be the solution matrix 
to the first order system corresponding to J 7 £ 7 [u]u = \iv written in the basis Yi(x) defined in (11), 
arguments similar to those above yield for small e 

( 0{\e\) 0(\e\) 

det(W(T;0,e)-W(0;0,e))= 0(\e\) V'{u-)T a + 0(\e\) V'{u-)T E + 0(\e\) 

V 0{\e\) 0{\e\) 

and hence D o (0,s) = 0(|e| 2 ) as claimed. □ 
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We are now in position to prove our main result of this section. By the above work, we can apply 
Lemma 3 to the equation D (n,s) = 0. The next theorem uses the Frcdholm alternative to discount 
possibilities (i) and (ii) from cmma 3, and establish the analyticity of the eigenvalues near \i = 0. 
Basically this amounts to checking that (generically) the null-space of the linearized operator has the 
same Jordan structure as the monodromy map at the origin. 

Theorem 4. For small e, the linear operator J 7 £ 7 [w] has three eigenvalues which bifurcate from \i = 
and are analytic in e. 

Proof. The idea of the proof is to systematically discount possibilities (i) and (ii) from Lemma 3, thus 
leaving only the third possibility. First, suppose case (i) holds. It follows from the Dunford calculus that 
we can expand the eigenvalues and eigenfunctions as 

f A*-e 1/2 ^i+e^+0(|e| 3/2 ) 

I « = /o+e 1/2 /i+£/ 2 +0(|e| 3/2 ) 

We will show the assumption that {T, M, P} a ,E,c 7^ implies v\ = 0, which yields the desired contradic- 
tion. Using the above expansions of v, (J, and J 7 £ 7 [</>] in terms of e, the leading order equation becomes 
Lafa = 0. Thus, / = &o0o + fri0i f° r some bo, b\ G C. Continuing, the Od^ 1 / 2 ) equation turns out to be 
Lof\ = vifo- Suppose vi ^ 0. By the Fredholm alternative, this equation is solvable in H if and only if 
&o0o + bi(f>i -L N(Lq). Clearly, <j>± _L N{L\) since 4>\ € Range(L ). Moreover, by Lemma 4 O ^ N(Ll) 1 - 
and hence we must have 6 = and, with out loss of generality, we take b\ = 1. It follows that f\ must 
satisfy the equation 

£o/l = ^101, 

i.e. /i = v\<$>2 + b 2 (po + &301 f° r some constants b 2 , 63 € C. 
Continuing in this fashion, the C(|e|) equation becomes 

£0/2 = ^1/1 + ^2/0 - Lxfo- 
By the Fredholm alternative, this is solvable if and only if 

(■00, "lfi + v 2 fo - Lifo) = and 
(ip2,vifi + v-ih - L i.fo) = 0. 

By above, /o is an odd function and since L\ preserves parity, the solvability condition implies we must re- 
quire (V>cb /1) = (02, fi) = 0. However, this is a contradiction since (ip 2 , <fa) = 2 {T, M} a ^{T, M, P} a ,E,c 
and hence it must be that v\ — as claimed. Thus, possibility (i) can not occur. 

Next, assume case (ii) of Lemma 3 holds. Then the Dunford calculus again implies the eigenvalues 
and eigenvectors can be expanded in a Puiseux series of the form 

f M = ^i£ 2/3 +^2£ 4/3 + 0(|e| 2 ), 

\ v = w + e 2 ^ Wl + e A ^w 2 + 0(\e\ 2 ). 

Our goal again is to prove the assumptions {T,M,P} a ,E, c 7^ and {T,M} a ^E 7^ imply u\ = 0. 
Substituting these expansions into J 1 C 1 [u]v = fiv as before, the leading order equation leads to wo — 
a o 0o + a i0i and the 0(|e| 2 / 3 ) equation implies a = 0. Without loss of generality, we assume a\ — 1, so 
that it follows that w\ = wi02 + a20o + as<t>i- The solvability condition at 0(|e| 4 / 3 ) implies that 

-Uj{ (02,02) = 0, 



which implies u)\ — as above. Thus, case (ii) of Lemma 3 can not occur leaving only case (Hi), which 
completes the proof. □ 
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4.3 Perturbation Analysis of spec( J 1 C 1 [u\) near (//, 7) = (0,0) 

We are now set to derive a modulational stability index in terms of the conserved quantities of the gKdV 
flow. By Theorem 4, it follows that the eigenvalues and eigenvectors are analytic in e, and hence admit 
a representation of the form 

J v = v + v 1 e + v 2 e 2 +C(|e| 3 ), 
\ n= A l£ + A 2£ 2 +0(| e | 3 ). 

At this point, it is tempting to use the functionals Pj := (tpj, ■) to compute the matrix action of the 
operator J 7 £ 7 [u] onto the corresponding spectral subspace associated with N g (Lo). This would convert 
the above eigenvalue problem for a fixed 7 to the problem of solving the polynomial equation 

det [Mo + eMi + e 2 M 2 + 0(e 3 ) - XP] = 0, 

at 0(e 2 ), where = {PiLkfa}ij and P = {Pifa}ij- Although this approach has been used to determine 
stability in the case where the underlying periodic waves are small (see [16] and [21]), this approach is 
flawed in the current case since, as shown below, the eigenvector v has a non-trivial projection onto 
Ng(Lo) of size O(e). Since we have no information about what such a projection would look like, it is 
unlikely that one can determine the nature of the spectrum near fj, = by computing the matrix action of 
the operator J 1 C 1 [u] on H for a general periodic solution of (1). Instead, we proceed below by developing 
a perturbation theory for such a degenerate eigenvalue problem based on the Fredholm alternative. 

Substituting the analytic representation of the eigenvector and eigenvalue into the equation J 7 £ 7 [u]v = 
[iv, the leading order equation implies vq G N(L n ), i.e. vo = c n fa + c\fa f° r some c , c\ G C. At C(|e|), 
we get the equation L vi = (Ai — L\)vo, which has corresponding solvability conditions 

= (ip ,L Vi) = Aic (tpo, <t>o) ~ Co (tpo, Li<j) ) - ci (-00, l>\§\) , and 

= (lp 2 ,L Vi) = -C (tp2,LKt)o) - Ci {lj)2,Li4>i) . 

It follows that we must require c = 0. Indeed, from the parity relation (ipi, Lk<f>j) = if i + j + k = 
mod(2), and the relations (ipo, (j) ) = {T, M} a ,E 7^ and (tpo,Li<fio) = Te ^ 0, we either have c = or 
all three eigenvalues bifurcating from \i = have the same leading order non-zero real part, which is not 
allowed by the Hamiltonian symmetries of the spectrum (recall spec(d x £[u]) is symmetric about the real 
and imaginary axis). With out loss of generality, we then set C\ = 1 and fix the normalization 

(iP u v) = (Vi,vo) = -\{T,M} aiE {T,M,P} a , EiC 
for all e. It follows that vo = <f>i and hence v\ satisfies the equation 

Lovi = (Ai - L^fa. 

Notice that L\Vq = — 2{T, M\ a ^u xxx does not belong to N g (L ), and hence the eigenfunction v has a 
non-trivial projection onto Ng(Lo) 1 - of size 0(e), as claimed above. 

We now define Lq 1 on R(Lq) with the requirement that R(Lq 1 ) is orthogonal to span{i})Q, This 
requirement ensures that Lq 1 f is well-defind and unique for all / € R(Lq). In particular, it allows us 
to compute the projection of L^ 1 f onto N(L a ) for each / € R(L ). In order to express the explicit 
dependence of V\ on Ai, we now write 

vi = L^ 1 (Ai - Li)4>x +c 2 (f>o +c 3 0i (22) 

for some Ci, C3 G C. The above normalization condition implies (ipi,vi) = 0, i.e. 

= (Vi, Lq \\! - Li)0i) + c 3 (Vi, fa) ■ 
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It follows C3 = by the definition of L 1 and the fact that (^1, <f>i) 7^ 0. 
Continuing, the C ( | e | 2 ) equation is 

io«2 = —Lwi - L 2 v + Xivi + X 2 v 

with corresponding solvability conditions 

= - {i/>o,LiVi) - (ipo,L 2 v ) + Ai (ipo, v-i) , and 
= - {ip2,LiVi} - (tl> 2 ,L 2 v ) + Ai {i[) 2 , vi) . 

Using the explicit dependence of vi on Ai and c 2 , it follows that we can express the above solvability 
conditions as 

Pi(Ai) + Pi(Ai)c 2 = Oand 
P 2 (Ai) - P (Ai)c2 - 0. 

As this is an over determined system of linear equations for c 2 , the consistency condition 

P(Ai) := P (Ai)Pi(Ai) + Pi(A 1 )P 2 (A 1 ) = 

must hold. In particular, this expresses Ai as a root of a cubic polynomial with real coefficients. Since e is 
purely imaginary, modulational stability follows if and only if -P(A) has three real roots, and hence it must 
be that A(/; u) is a positive multiple of the discriminant of the cubic polynomial P{\). Notice that one can 
explicitly calculate P(A) for a general non-linearity using just the definitions of the <j>j and ipj, except for 
the inner products (Vo, L\Lq 1 Li0i) and (ip 2 , LiL^ 1 L\(j)\) . The first of these can be calculated regardless 
of the nonlinearity, but we must restrict to power-law nonlinearity for the computation of the second (see 
appendix). It follows that we can explicitly write down the compatibility condition P(Ai) = only in 
terms of the underlying periodic solution u and terms built up out of the generalized null spaces of Lq 
and Lq acting on L 2 (R/TZ). Since the roots of this polynomial determine the structure of spec(J 7 £ 7 [u]) 
in a neighborhood of the origin, we have proven the following theorem. 

Theorem 5. The periodic solution u = u(x; a, E, c) of (1) is spectrally unstable in a neighborhood of the 
origin if and only if the discriminant A(a, E, c) of the real cubic polynomial P(A) is positive. Recall that 
the discriminant of a cubic P(A) = aX 3 + b\ 2 + c\ + d is given by A = & 2 c 2 — 4ac 3 — 46 3 rf — 27a 2 d 2 + 18abcd 

Remark 6. The above result gives a second characterization of the modulational stability of periodic 
solutions to the generalized Korteveg-DeVries equation with power law nonlinearity since it is expressed 
entirely in terms of T, M, P, H and their derivatives, which in turn can be written as functions of a, E, c 
via integral type formulae. (These are hyperelliptic integrals in the case that p is rational). The formulae 
remain, however, somewhat daunting. Since this detects the same instability that the Evan's function 
based criterion does this quantity must have the same sign as the discriminant derived in that section, 
although we have not shown this explicitly. 

5 Concluding Remarks 
5.1 Discussion 

We'd like to consider a concrete example to illustrate our results. We have chosen to consider the power 
law gKdV with p = 5(andc =1). In this case the solitary wave is unstable and hence (by Gardner's 
result, which we have checked in this case using our methods) periodic waves of sufficiently long period 
are also unstable. Haragus and Kapitula[21] have done some very nice experiments on this case using 
the SpectrUW[12, 13] package, which they have been kind enough to share with us. For clarity we have 
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Figure 3: Cartoon of the spectrum of the linearization of gKdV about a periodic traveling wave for 
p = 5, a = 0, E = and three different values of the period (ordered by increasing period). 

illustrations representing the spectra they computed numerically rather than reproducing their figures - 
see Figure 3. 

The first graph in Figure 3 depicts the spectrum for small amplitude periodic waves (this is the solution 
branch inside the separatrix). The modulational instability index A < indicating a modulational 
instability while the orientation index {T,M,P} aEc > 0. The latter indicates that the number of 
eigenvalues on the real axis away from the origin is even. In this case there are none. The spectrum near 
the axis looks like a union of three straight lines, as predicted by the fact that the normal form of the 
periodic Evan's function is a homogeneous polynomial of degree three. Globally the spectrum looks like 
the union of the imaginary axis with a figure eight shaped curve. 

As the period increases one sees spectra which resemble the second figure, where there is a modula- 
tional instability together with a pair of eigenvalues along the real axis. In this case we are still in the 
case A < 0, indicating a modulational instability and {T, M, P} a ,E,c > indicating an even number of 
eigenvalues along the positive real axis. The fact that these two very different spectral pictures have the 
same orientation and modulational instability indices shows that these quantities are not enough to say 
qualitatively what the spectral picture looks like, even in this very simple problem with only one free 
parameter (the period). 

As the period increases still further one sees spectral pictures which resemble the third picture. 
As in the previous figure there is an oo shaped curve of spectrum connected to the origin indicating a 
modulational instability (A < 0) as well as two loops of spectrum intersecting the real axis and supported 
away from the origin. These loops are those predicted by Gardner in his paper arising from the discrete 
eigenvalues of the solitary wave problem. As the period increases and the periodic solution approaches 
the solitary wave the circle collapses to a point and the oo curve collapses to the origin. The size of both 
of these features is exponentially small in the period. In the paper of Kapitula and Haragus the oo curve 
is not visible at the scale of the graph, but it is visible in numerics they performed for smaller values of 
period. 

Since there is an odd number of eigenvalues on the real axis in this case (one periodic, two antiperiodic) 
the orientation index must now be negative: {T, M, P} a ,E,c < 0. The general mechanism by which this 
must occur is clear: a periodic eigenvalue moves down the real axis, collides with the origin (changing the 
Jordan structure of the null-space of the linearized operator, which is again signalled by the vanishing of 
{T, M, P} a ,E,c <) and moves off along the real axis. However the exact way in which this occurs is not 
quite clear. 
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5.2 Open Problems and Concluding Remarks 

We have considered the stability of periodic traveling waves solutions to the generalized Kortcvcg-DeVries 
equation to perturbations of arbitrary wavelength. We introduce two indices related to the stability of the 
period waves. The first, which is given by the Jacobian of the map between the constants of integration 
of the traveling wave ordinary differential equation and the conserved quantities of the partial differential 
equation, serves to count (modulo 2) the number of periodic eigenvalues along the real axis. This is, 
in some sense, a natural generalization of the analogous calculation for the solitary wave solutions, and 
reduces to this is the solitary wave limit. The second, which arises as the discriminant of a cubic which 
governs the normal form of the linearized operator in a neighborhood of the origin, can also be expressed in 
terms of the conserved quantities of the partial differential equation and their derivatives with respect to 
the constants of integration of the ordinary differential equation. This discriminant detects modulational 
instabilities: bands of spectrum off of the imaginary axis which are connected to the origin. As we have 
emphasized throughout this calculation can be considered to be a rigorous Whitham theory calculation. 

This calculation hinges on the fact that the underlying ordinary differential equation has sufficient 
first integrals. As such it is doubtless related to the multi-symplectic formalism of Bridges [9] . As there are 
a number of other equations for which the traveling wave ODE has a integrablc Hamiltonian formulation 
(Bcnjamin-Bona-Mahoney, NLS, etc) one should be able to carry out the analogous calculation in those 
cases. The additional structure provided by the scaling invariance is also extremely helpful, as this allows 
one to simplify many of the calculations but, at least in the Evan's function approach, it has not been 
necessary. 

We are somewhat puzzled by the fact that the Evan's function based calculation gives a substantially 
simpler criteria for the existence of a modulational instability than one based on a direct analysis of the 
linearized operator. It must be true that the two discriminants we've derived always have the same sign, as 
the predict the same phenomenon, but we have been unable to see this directly from the formulae. Often 
when apparently unconnected quantities share a sign this sign has a topological or geometric interpretation 
(for example as a Krein signature), so this may well be the case here. Such an interpretation would be 
very interesting. 
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6 Appendices 
6.1 Identities 

In this section we derive a number of useful identities which will allow us to relate various Jacobians 
which arise in the calculation. Wc define the conserved quantities: 



r 

T=(l)=2 

J U 

M = {u)=2 
P = (u 2 ) = 2 



u+ 



du 



^2E + 2au - 


- -pp2uP+ 2 + cu 2 


+ 


udu 


^2E + 2au 


- ^u p+2 + cu 2 


*+ 


u 2 du 



H = 



ul U p+ 2 



^2E + 2au - ^uP+ 2 +~^ 2 

u + (E + au-^ + cu 2 ) 



P + 2 



2E + 2au - -pp2 uP+2 + cu 2 



The classical action K = § u x du = u 2 dx provides a useful generating function, and is given by 



K = 2 



rUJ, 
J U- 



2E + 2au - 



P + 2 



u p+2 _|_ cu 2 d u . 



This integral has the advantage that it is regular at the endpoints and can thus be differentiated in the 
form presented. It is obvious that the derivatives are given by 



dK 

da 
dK 

~dE 

dK _ P 

~dc~^2 



= M 
= T 



(23) 



Note that these relations force certain relations among various 2x2 Jacobians. For example we have 



{M,P} a , E 
Similarly we have the relations 



M a 


Pa 






M E 


Pe 




K aE 2K cE 



K aa KaE 
K ac K cE 



= -2{T,M} a ,, 



{T,M} E . c = --{T,P} a . E 
{T,P} a , c = {M,P} E , c . 

There is another identity relating the gradients of T and the conserved quantities M, P, H which is useful. 
Begin by noting that 



2 \p + 2 

1 -k-(^) n 

2 \p+2/ 



= ET + aM + -P 
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c 
2 









V i 


J 




V o / 



where the former quantity is the Lagrangian of the traveling wave differential equation. Adding these 
together, taking partial derivatives with respect to (a, E, c), and using the relations (23) shows that 



T a M a P a H a 
T E M E P e H e 
T c M c P r H r 



so that the gradients of any of T, M, P, H can be expressed in terms of an explicit linear combination of 
the other three. As noted in the text the theory is most developed in terms of the first three quantities, 
as these arise most naturally, but is stated in terms of the last three, as these have the most natural 
physical interpretation. 

There is another set of Jacobian identities which are useful. Differentiating (2) with respect to E and 
subtracting ue times (5) gives the identity 



Similarly we have the identities 



u x u xE - u xx u E = {u, u x } x . E = 1. 



{u,u x } x . a = u 
{li, u x ^ x c —u . 



There are a number of other identities of this sort which can be derived in an analogous fashion. 



6.2 Analysis of N(C[u]) 

In this appendix, we give a detailed analysis of the null space N(£[u]). As above, we assume u is a 
solution of (1) of period T and define H = L 2 (R/TZ;C). Notice that from Proposition 2 we know there 
that the functions u x , ue both satisfy the differential equation = when boundary conditions are 

ignored. However, u E is not in general T-periodic due to the variation in the period with respect to E. 
This gives an indication that u x is the only T-periodic solution of £[u]v = unless one has Te = 0. This 
observation leads us to the main result of this appendix. 

Theorem 6. Considered as an operator on H, N(£[u]) = span{u x }. In particular, up to constant 
multiples, the equation £[u]v = has only one solution in H. 

Proof. Define y\(x) = (^jff^ ue{x) and 2/2(2) = — V'(u-)~ 1 u x and note that 

yi (0) = 1 y 2 (0) = 
2/i(0) = y' 2 (0) = l. 

Moreover, £[u]yj = for j = 1, 2. In the y\, y 2 basis, an easy calculation then proves the monodromy is 
expressed as 




and hence is a band edge of spec(£[w]). It follows that there exists a second periodic element of 
N(£[u]) if and only if Te = 0. The proof that Tg > is a result of the following theorem by Schaaf, 
which completes the proof. □ 



6 APPENDICES 



30 



Lemma 5. Assume G is aC 3 function on (0, oo) and that G vanishes only at one point Xq with G'(xq) > 
0. Define 

A = {x e M : x < x and G(s) < for all s e (x, x )} 
and suppose for each a € A there exists a periodic solution x(t) > of the equation 

x"(t) + G(x(t)) =0 

with initial data x(0) = a, x 1 (0) = 0. Let P(a) denote the period of this solution If G satisfies the two 
conditions 

G'{x) > 0, x g A =4> 5G"(x) 2 - 3G'(x)G"'(x) > 0; 
G'(x) = 0, x e A => G(x)G"(x) < 

then P is differ entiable on A and P'(a) > 0. 

In order to apply the above result in our case, define G(x) = V'(x;a) and assume G(x a ) = 0, 
G'(x a ) > 0, and c = 1. Define A as above and notice that G € C 3 (0,oo). For all x such that G'(x) > 0, 
we have x p > and hence 

(5G" 2 -3G'G"')(x) = P (p+l)xP- 2 [(p+l)(2p + 3)xP + S(p-l)} 
> hp 2 {p + l)x p - 2 > 

for such x. Also, if G"(a;) = 0, then x p — ^j-j- and hence at such points, 



G{x)G"{x) =p 

given that 



p+l 



1 -a(p+ 1)^- 1 



< 



aa;^ 1 = a I I 



p-i 

1 1 



,P+1J (P+l) 2 P+l 

for all x e A. Hence, it follows for such a that Tg > for all periodic waves bounded by a homoclinic 
orbit in phase space. 

Next, as pointed out in the text, the fact that the monodromy matrix at the origin is the identity 
plus a rank one perturbation implies there is a linear combination of u a and ue which is periodic. This 
combined with the above lemma immediately implies the following. 

Corollary 3. Considered as operators on H, we have N(d x C[u]) = span{(j) , 4>i} and N(C[u]d x ) — 
span{l, u}. 

We end our discussion with the following interesting remark. By our above work we see that the 
requirement Te 7^ is sufficient of the origin to not be a double point for the sepctrum of the Hill- 
operator C[u]. There is a geometric quantity which detects this same information known as the Krein 
signature. For the operator £[u], the Krein signature at the origin is easily shown to be tr(m'(0)), where 
m(0) is defined as above. To see this, notice the characteristic polynomial for m(p) can be expressed as 

det[m(p) - XI] = A 2 — tr(m(ju))A + 1 

which roots 

tr(m(/i)) ± v /tr(m(/x)) 2 -4 

A± = 2 ' 

Thus, the solutions to the equation tr(m(/i)) = ±2 correspond to the periodic eigenvalues of the operator 
£[u] and, moreover, p is a double point of the periodic spectrum if and only if tr(m'(/i)) = 0. From 
this discussion, it follows that tr(m'(0)) must be a non-zero multiple of Tg, a fact which is proven in the 
following proposition. 
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Proposition 5. For the operator C[u], one has sgn (tr(m'(0)) = sgniTE)- As a result, is never a 
double point o/spcc(£[u]) under the assumption T E ^ 0. 

Proof. This proof is essentially given in Magnus and Winkler[29]. Using variation of parameters, one can 
express j^y.j and -j^y'j in terms of the yj and y'j. Using the facts that det(M(0)) = 1 and tr(m(0)) = 2, 
a bit of algebra eventually yields the expression 

tr(m'(0)) = sgn(y[(T)) ^ U\y' x {T)\y 2 + S gn(y[(t))^l=^l yi \ dx 



It follows that 

sgn(T(m'(0))) = sgn(y[(T)) 



( du_ 



-l 



= sgn(T E ) 

as claimed. □ 
6.3 Negativity of M a 

In this appendix we show that M a (a,0) < for a sufficiently small and c = 1. Note that M(a,0) can 
(after some rescaling) be expressed in the form 



r r(a) 

\/a + U — UP+ 2 



M(a,Q)= f v = ^ du 

Jo 



where r(a) is the smallest positive root of a + u — u' p+2 = 0. This is a smooth function of a for a small 
enough and satisfies 

r(a) = l + ^ + 0(a 2 ) 

The main idea is to rescale the above so that the integral is over a fixed domain and show that the 
integrand is a decreasing function of a on the new domain. Rescaling gives 



M(a,0)= / v du 

The quantity ^ + ^ - r^u^ 2 satisfies 

4 + ^ - r-P- V+ 2 = u - U p+ 2 + o(l - —v - ^-v?+ 2 ) 

r A r Z p + 1 p + 1 

The second term is clearly positive on the open interval (0,1), and thus the rescaled integrand is a 
decreasing function of a, and thus M a < for a small enough. 

6.4 Evaluation of Virial-Type Identities 

In this section, we evaluate the two virial type inner products arising in the perturbation analysis of 
section 4.3. One of these is calculable for an arbitrary nonlinearity, while for the other we must restrict 
to the case of power-law nonlincarities. We proceed with the more general one first. 



Lemma 6. {^L^L^) = -T{T,K} a . E . 
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Proof. Define an operator £ : {g € L 2 (R/TZ) : (g) ^ 0} -» L 2 (R/TZ) by 

ag)=x ~v)L g{s)ds - 

Then a straight forward computation shows that Lq£(0o) = f'( u ) — c + {t(m)T e ' ^ f° uows that 

(Vciiio^i^i) = 2{T,M} 0>B ((/'(«) -cJ.Lo^xxx) 

= 2{T,M} a . E ^LlaM,L^u x: 

= T((f) ,u xx } 

= —T{T, K} a e 

as claimed. □ 

While the above expression holds for an arbitrary nonlinearity, we have found a closed form expres- 
sion of (ip2, LiL^ 1 Lo4>i) in the case of power non-linearities. From the evaluation of the modulational 
instability index via Evans function techniques, it should be that this inner product is calculable in the 
general case as well, although we have yet to be able to do this. 

Lemma 7. In the case of a power nonlinearity f(x) = x p+1 , we have 
{^LxL^L^) = -T{T,M} E ,c{T,K} a , E 

+ 2 —{T, M} a e (M{T, K} a , E - 2{T, M} a , E K) 
P 

+ 2c{T,M} aiE {T,M,K} a!E!C . 

Proof. Notice that in the case of power-law nonlinearity, one has 

(^LiL^L^) = -T{T,M} E . c {T,K} a .E 

+ {T, M} a . E ((2 - p) (u p+1 ,Lo 1 L^ 1 ) - 2c (u, L^L^)) , 

and hence we must evaluate (u, L$ and (u p+1 , L^ 1 Li4>i) . First, from the definition of v\ in 

equation (22) it follows that 

(tp2,vi) = Ai (V>2,io Vi) - (^2,^0 lL i<t>i) 

= ~\i{T, M} a , £ {T, M, P} a , E ,c ~ {T, M} a E (u, L„ 

Moreover, using the fact that ip2 — L' q i/ji gives 
(ip2,vi) = (ipi, (Ai - 

= ~\i{T, M} Qi£ {T, M, P} atE ,c - 2{T, M} a , E (<h, u xx ) 

= ~\i{T, M} a , £ {T, M, P} a , E ,c + {T, M} a>E {T, M, K} a , E , c 

and hence (u, L~ 1 L 1 cf) 1 ) = -{T, M, K} a . E . c - 

Next, let the functional £ be as in Lemma (6) and notice that 

LU(u) = f(u)-c-^(pu p +'+a). 
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It follows that 

-^<^+ 1 ,L-%0 1 ) = (4£( U ) _(/'(„) _ c ), L -i Ll ^ 

= (£(u), Li^i) + 2{T, M} a , £ ((/'(«) - c),L^u xxx ) 

2T{T,M} a , E T 
= K - l{L,K\ a . E 

which completes the proof. □ 
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